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The conditions under which ferromagnetism may occur in transition metal oxides with partially 
filled t2g shells such as vanadium-based perovskites are studied using a combination of density 
functional and single-site dynamical mean field methods. For reasonable values of the correla- 
tion strength, rotations of the VOe octahedra play an important role in enabling ferromagnetism, 
with ferromagnetism typically occurring for rotations larger than a nonzero critical value. Fer- 
romagnetism is suppressed near the Mott-insulating phase but the phase boundary is otherwise 
only weakly dependent on carrier concentration. Design rules are suggested for new oxide systems 
exhibiting ferromagnetism. 
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I. INTRODUCTION 

Transition metal oxide^ are of great interest in con- 
densed matter physics because they exhibit a rich va- 
riety of exotic phenomena which remain incompletely 
understood. While "late" transition metal oxides (in- 
volving Cu or Ni) have been very extensively stud- 
ied due to their connection to high- Tc superconductiv- 
ity and Mn-based compounds have attracted attention 
for their colossal magnetoresistance, the "early" transi- 
tion metal oxides such as vanadium oxides have, despite 
some important studies,'^ received less attention in re- 
cent years. However, following the pioneering work of 
Ohtomo and Hwang,!^ early transition metal oxides are 
increasingly used as components of atomic-scale oxide- 
based heterostructures.'^ One goal of heterostructure 
research is to design materials exhibiting phases not ob- 
served in bulk.-^ An essential step towards realizing this 
goal is obtaining a clear understanding of the relation- 
ship between physical structure and observed electronic 
phenomena. 

In this paper we investigate the relationship between 
lattice structure, correlation strength and electronic 
properties in the context of ferromagnetism in early tran- 
sition metal oxides. Ferromagnetism is a correlated elec- 
tronic property which is both technologically important 
and (because only a spin symmetry and not translation 
or gauge symmetry is broken) more straightforward from 
the theoretical and computational points of view than 
other phases such as antiferromagnetism or supercon- 
ductivity. We choose the early transition metal oxides 
in part because of the intriguing recent reporiP of ferro- 
magnetism in superlattice systems involving LaVOa and 
SrVOa . The report is of interest because ferromagnetism 
is reported for the superlattice even though no ferromag- 
netism is observed in bulk solid solutions of the form 
Lai_a:Sra.V03. One possible explanation is that the su- 
perlattice enables a different crystal structure, more fa- 
vorable to ferromagnetism than that found in the ob- 
served bulk structures. Understanding whether this ex- 
planation is viable, and more generally being able to de- 



sign superlattices with desired magnetic properties, re- 
quires deeper insight into the conditions 

The conditions under which ferromagnetism may oc- 
cur is a question of long-standing theoretical interest .ISHni 
The development of dynamical mean field theorjff^l has 
opened a new avenue of research, but apart from some 
pioneering investigations of the Curie temperatures of Fe 
and Np31 the studies have mainly been based on model 
systems. VoUhardt, Ulmke, and collaborators have stud- 
ied the single-band Hubbard model, finding that in this 
model ferromagnetism occurs at generic carrier concen- 
trations only when there is a lar ge den sity of states peak 
at or near the lower band edge^^^^ For a fixed value 
of the Hubbard U, the Curie temperature Tc was found 
to depend sensitively on the peak position, becoming un- 
measurably small as the density of states peak was moved 
a small distance away from the lower band edge.l^^ How- 
ever, many materials of physical and technological inter- 
est involve transition metals with partially filled degen- 
erate (or nearly degenerate) d levels, where the Hund's 
interaction may play an important role. While the im- 
portance of the Hund's interaction in partially filled d 
levels has been appreciated for decades, the issue has 
been systematically studied only in the case of the Bethe 
lattice^^Kl^ in which the density of states has a simple 
semicircular structure. In this situation large values of 
the Hund's coupling and correlation strength (the Hub- 
bard U) are required for ferromagnetism. 

In this paper, we study the conditions under 
which ferromagnetism may occur in systems related to 
Lai-ajSr^^VOa (LSVO). In this family of materials the key 
orbitals are vanadium d levels. These are partly filled (in 
Lai_a;Sr2:V03 the mean number of d electrons is 2 — a;) 
and the electrons are subject to a strong Hund's coupling, 
favoring magnetism. The end-member LaVOs is classi- 
fied as a Mott insulator.'^ It has structural and antiferro- 
magnetic transitio ns at around 140K but is insulating to 
high temperatures.'2212ll Our investigations are based on 
density functional theory (DFT) calculations and include 
in particular the effects of GdFeOa-type rotations of the 
BOg structural motif of the ideal perovskite structure. 
We treat the many body physics using the single-site dy- 
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FIG. 1: Three possible magnetic phase diagrams for bulk 
solid solution Lai_a:Sra:V03 as a function of carrier density 
{x axis) and tilt angle {y axis). The phase boundaries are the 
black curves, (a) Weak correlations: Slater antiferromagnet 
at a; = 0. (b) Intermediate correlations, (c) Strong correla- 
tions: Mott insulator at a; = 0. The notations are: FM for 
ferromagnetic, FM for paramagnetic, AFM for antiferromag- 
netic order, and I for insulating state. 



LaVOa strongly increases the tendency towards ferro- 
magnetism. An intermediate situation may also occur. 
The three panels of Fig. [T] present schematic representa- 
tions of these three cases. 

We remark that in the literature, the many-body 
physics properties of transition metal oxides have been 
modeled in two ways. One, which we refer to as the "d 
only model," treats the frontier orbitals (in the case of the 
vanadium-based materials, antibonding bands derived 
mainly from V t2g-symmetry d orbitals admixed with 
oxygen states) as a multiorbital Hubbard model. This 
approach is generally accepteciD as the relevant descrip- 
tion of the early transition metal perovskites such as the 
systems we study here. A second approach is based on 
atomic-like d orbitals which are subject to on-site corre- 
lations and are hybridized with the full p manifold. This 
approach is used in the DFT+U and DFT+DMFllSami 
approaches. It is generally believed to be essential to 
a correct description of the late (e.g. Cu-based or Ni- 
based) transition metal oxide^ but has been less widely 
used in the study of early transition metal oxides. We 
consider both approaches in this paper, finding similar 
qualitative conclusions but significant quantitative dif- 
ferences which arise from the richer physics of the "pd 
model" situation. 

The paper has following structure. The model and 
methods are described in Sec.jlll Section ITTTl derives the 
magnetic phase diagram of LaVOs as a function of carrier 
concentration, tilt angle and correlation strength using 
the widely used d only model. In this analysis, most at- 
tention is given to the strongly correlated case presented 
in Fig. [lb. Section |IV| analyzes the pd model. The fi- 
nal Sec. |V| is the summary and conclusion. Appendices 
present details of some of the results. 



namical mean field approximation'-'^ - which is wi dely u sed 
in real-materials many-body physics calcul ation j^^^ and 
in particular in studies of fcrromagnetism.'SIlll 

Our main result is an understanding of the dependence 
of ferromagnetism on carrier concentration, octahedral 
rotation and correlation strength. At non-integer elec- 
tron concentration we find that increasing the octahe- 
dral rotation amplitude increases the tendency to fer- 
romagnetism. The physics underlying this result is a 
density of states effect related t o that previously noted 
in the one-band Hubbard model.'^^'^ll We also find that 
at fixed rotation amplitude the carrier concentration de- 
pendence of ferromagnetism depends on the correlation 
strength. For weak to moderate correlations (so that at 
the carrier concentration n = 2 corresponding to LaVOa 
the paramagnetic phase is metallic in contradiction to 
experiment! ^ l ^° l ^^ l) decreasing the carrier concentration 
below n = 2 weakly decreases the tendency to ferromag- 
netism. For stronger correlations (so that at the carrier 
concentration n — 2 corresponding to LaVOa, the in- 
sulating phase is a Mott insulator), doping away from 



II. MODEL AND METHODS 
A. Model 

We study systems derived from SrVOa and LaVOa. 
Bulk SrVOa crystallizes in the ideal cubic ABO3 per- 
ovskite structure.'2£l The crystal structure of LaVOa may 
be thought of as a cubic perovskite with an additional 
GdFeOa-type rotational distortion leading to a tilted 
structure with symmetry group Pnma-^ Fig. [2] shows 
the lattice structure of LaVOa. The basic structural mo- 
tif of the perovskite structure is the oxygen octahedron. 
In the GdFeOa rotated structure there are four inequiva- 
lent octahedra characterized by different directions of the 
principal axes and by slight (~ 5%) distortions of the V- 
O bond length. We have found (not shown) that varying 
the relative magnitudes of the V-0 bond lengths over 
the physical range causes only small changes to the non- 
interacting DOS or to the DMFT solution. We therefore 
set all V-0 distances to d = 1.95A and focus on the 
effect of rotation by studying a range of 6 and 

The octahedral rotations in the GdFeOa-distorted 
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online): A atoms (here La or Sr), small spheres (red online): 
O atoms, intermediate size spheres at the centers of the octa- 
hedra (blue online) : B-site atoms (here, V). 



Pnma systems can be characterized by two angles, 9 
and 0,'^ with corresponding rotation axes fig, fi^ and 
wavevectors Qo, characterizing the changes in rota- 
tion axis from cell to cell. We choose coordinates such 
that the rotation axis hg is [110] while is [001]. The 
corresponding wavevectors are Qg = (tt, tt, tt) and = 
(7r,7r,0). For bulk LaVOg, 6 = 11.5° and = 8.8°.^° As 
La is progressively replaced by Sr in bulk solid solution 
LSVO, 9 and (p go to 0. Appendix A provides an esti- 
mate of the x-dependence of 9 and <p hi the actual solid 
solution. 

We now turn to the electronic structure. The large 
energy splitting between the transition metal d bands and 
oxygen p bands characteristic of early transition metal 
oxides such as La/SrVOs is generally believecPto justify 
a downfolding to a model containing only the frontier 
bands, which in the present case are composed mainly 
of V-derived t2g symmetry d states with only a small 
admixture of oxygen p states. We refer to this as the 
d only model. In a subsequent section we present an 
examination of a more general "charge transfer" model 
in which the full p d complex is considered. 

The Hamiltonian of the d only model is 

H = Hkin + Honsite, (1) 

where Hkin describes the electron hopping between sites 
and Honsite describes the d d interactions, which we take 
to be on-site. 



The kinetic Hamiltonian has the quadratic form 

Hk^n = ^6'fnd(k)cLaCk/3,T, (2) 

k,Q,/i^,cr 

where cj^^,^ and Ck/So- are electron creation and annihila- 
tion operators in reciprocal space with wavevector k. a 
and P are orbital indices, and a is the spin index. 

For the cubic structure, the HbandO<) describing the 
^2g-derived antibonding bands is easily constructed using 
phenomenological tight binding methods.l^For example, 
for calculations with only antibonding bands composed 
mainly of t2g orbitals, HbaiidO^) is almost diagonal and 
well approximated by a two-dimensional energy disper- 
sion with nearest neighbor hopping t and next-nearest 
neighbor hopping t' . The explicit form for the band aris- 
ing from the xy-symmetry orbital is 

e(k)j;j, = — 2t(cos kx + cos ky) — At' cos k^ cos ky. (3) 

and for the xz / yz-derived bands one relabels the momen- 
tum indices appropriately. 

For the tilted structure, the form of HtandO^) is more 
difficult to express in the simple tight binding language 
used for the cubic perovskite structure. First, the natural 
basis for defining local d orbitals is aligned with the axes 
of the local VOg octahedron and so varies from site to 
site making the Slater-Kanamori procedure much more 
involved. Second, the V-O-V bond is buckled, decreas- 
ing the hopping t both in absolute terms and relative to 
t' , as well as opening many other channels of hopping 
from one V-site to another. For tilted structures we con- 
struct HbandO^) numerically from density functional the- 
ory (DFT) band structure calculations using maximally- 
localized Wannier function (MLWF) techniques.'^ We 
have verified that when applied to the cubic structure 
the method reproduces the phenomenological tight bind- 
ing model discussed above. 

We found that in the resulting MLWF band structures 
(not shown), the bands derived from and oxygen p 
orbitals agree very well with the DFT results, but the 
eg antibonding bands do not because they are entangled 
with many other bands. However, the eg states derived 
from these bands are well above the Fermi level and, in an 
appropriate local basis aligned with the symmetry axes 
of a VOg octahedron, do not mix with t2g bands. They 
therefore do not affect the final results significantly. 

The on-site interaction term in the Hamiltonian, 
Honsite^ assumcd to take the standard Slater-Kanamori 

Honsite — U 'Y^nia-fUiai + Ui ^ nia^njpi + 
+ U2 ^ niacrnii3a + 
+ Jex 'PLt'^mi'^f3l^^o.l + 
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where U is the on-site intra-orbital interaction, U2 and 
Ui are on-site inter-orbital interactions for electrons with 
the same and different spin respectively and J^x and Jpair 
are the exchange and pair hopping parts of the interac- 
tion. For a transition metal ion in free space the rela- 



tions between these coefficients are: Je 



J, 



pair 



= J, 



Ui = U - 2J, U2 = Ui - J = U - 3 JP It is generally 
believed that the relations also hold for the t2g shell - and 
we make this assumption here. We consider several val- 
ues of U and J but focus most attention on the values 
U = 6eV ~ 22t and J = leV. These values of U and J 
are similar to but slightly larger than those used in re- 
cent paper^^UMl ^nd are chosen to reproduce the crucial 
feature of the phase diagram, which is that LaVOa is a 
Mott insulator. 



C. Determining the magnetic phase boundary 

In order to detect magnetic order, one can allow for 
broken symmetry states in the DMFT procedure and 
lower the temperature until an ordered state is reached. 
However, ordering temperatures are typically so low rel- 
ative to the basic energy scales in the problem that our 
simulations become prohibitively expensive. We there- 
fore compute the susceptibility which is found to have a 
Curie- Weiss form ~ T — T^. Extrapolation of 
to yields an estimate of the Curie temperature Tc- We 
interpret positive Tc as evidence of magnetism. To com- 
pute x~^, we add a uniform field H-a to our Hamiltonian, 
compute the magnetization m{H), verify that m is linear 

in H and then define = — ■ 



B. Dynamical mean- field theory 

To treat the on-site interaction terms we use single- 
site dynamical mean field theory (DMFT), [i^! which ob- 
tains the electron self-energy E(a;) (assumed to be site- 
local but dependent on spin and orbital indices) from 
the solution of a quantum impurity model. We use the 
hybridization-expansion CTQMC solveJ^ to solve the 
quantum impurity model. In the present case, the im- 
purity model is a three-orbital model specified by the U 
and J interactions and by a hybridization function fixed 
by the self-consistency condition 



S(w) 



(5) 



where {}ii denotes the on-site Green's function at vana- 
dium site i. 

In applying the procedure to the GdFeOs-distorted 
structure, it is important to formulate the impurity prob- 
lem corresponding to a given V site in a local basis 
aligned to the V-0 bond directions of the VOg octahe- 
dron at this site. The self energy derived from solving 
the impurity problem for each octahedron is then ro- 
tated back to the global basis to define the lattice Green's 
function used in the self-consistency condition. There 
are three advantages. First, the hybridization function 
is nearly diagonal in the local basis, which greatly re- 
duces the sign problem in the quantum Monte Carlo so- 
lution of the impurity model. Second, the t2g and Cg 
bands are clearly distinguished and well separated in the 
local basis, enabling removal of the Cg bands from the 
many-body problem (they can be neglected or treated 
in mean field) so the dynamical impurity model involves 
three rather than five bands. Third, the self-consistency 
conditions corresponding to the four inequivalent V-sites 
of Pnma LaVOa are related by rotations so only one im- 
purity model needs to be solved. 
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FIG. 3: (Color online) Inverse magnetic susceptibility = 
— - — ^ plotted vs. temperature T obtained from single-site 

DMFT solution to Eqs. and 3| with Hkin taken to have 
the form appropriate to a simple cubic lattice with nearest 
neighbor hopping t and second neighbor hopping t' = — 0.3t. 
Carrier density is fixed to n = 1.5 and applied magnetic field 
H = 0.005 (circles) and H = 0.01 (diamonds). The Curie 
temperature Tc = 0.234t is estimated by linear extrapolation 
of x~^ in the region where m is linear in H. Inset: Expanded 
view of the near-Tc region, together with the magnetization 
squared (m^) (left y axis) calculated by allowing for broken 
symmetry DMFT solution at _ff = which also shows that 
Tc « 0.222f. The parameters are t = 0.5, t' = -0.3t,U = 

16 = 32t and J =—. 

6 

Fig.[3]demonstrates the method on a model with a high 

Curie temperature, where the magnetic state can easily 

be constructed. This model is defined by Eqs. ([I] [s]) 

with Hkin taking the same form as in the simple cubic 

t' 

lattice but with t' chosen as — = —0.3, so that the sign 
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oi t /t is opposite to the sign implied by band structure 
calculations on the actual materials. For this unphysical 
sign of t /t the model has a ferromagnetic ground state 
with a high Curie temperature. We set carrier density 
n = 1.5 and calculate Tc in two ways: by lowering the 
temperature until ferromagnetic order is observed, which 
is at Tc w 0.222t; or by measuring x~^(T) for several 
values of T above Tc and linearly extrapolating to — 
0. The extrapolation shows Tc — 0.234i. The two values 
are very close. We conclude that extrapolating Tc from 
X~^(T) is a reliable way to determine whether the model 
exhibits ferromagnetic order. 



III. MAGNETIC PHASE DIAGRAM 

In this section, we study the magnetic phase dia- 
gram of model systems derived from the calculated band 
structure of LaV03 but with variable amplitude for the 
GdFeOs distortion and different values for the carrier 
concentration. 
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FIG. 4: Density of states for bulk LaVOa (d carrier density 
n — 2) calculated from DFT with different values of tilt angle 
9 and rotation angle 4> in which 6 = (f) together with the 
DOS of realistic LaVOs structure {9 = 11.5°, = 8.8°). A 
van Hove peak at the Fermi level develops as 9 and </> increase, 
and is well-formed a,t 9 = (j) = 14° and above. The bandwidth 
W decreases as the 9 and 4> increases. The Fermi energy is at 
0. 

Fig. |4] shows the evolution of the non- interacting den- 
sity of states with increasing amplitude of the GdFeOa 
distortion. If the structure is cubic, the energy disper- 
sion in the t2g manifold is almost two dimensional. The 
sign and magnitude of t' place the van Hove peaks of 
the non-interacting DOS at high energy (see top panel 
(6* = (/) = 0) of Fig. |4]). As the distortion ampli- 
tude is increased the lowest DOS peak shifts to lower 
energy and the bandwidth decreases. Previous litera- 
ture suggests a connection between peak position and 
ferromagnetismfSmi while standard Stoner theorj^ sug- 
gests that bandwidth may also be important. 



4.5 



s 
o 



3.5 



2.5 




0.05 



- W 

0.1 0.15 
temperature (eV) 



X'forW=2.4eV 
X"'forW=1.6eV 
R„forW=2.4eV 
R„forW=1.6eV 



0.25 



0.2 



0.15" 



0.1 



a* 

X 



0.05 



0.2 



Jo 

0.25 



FIG. 5: (Color online) Inverse magnetic susceptibility 
(closed symbols, black on-line; right axis) and Wilson ra- 
tio Rw (open symbols, red on-line, left axis) calculated for 
Lai_2,Sr3;V03 cubic structure as functions of temperature at 
carrier density n = 1.5 and on-site interactions U = 6eV, 
J = leV using two bandwidths, W = 2.4eV (circles and solid 
lines) and W = 1.6eV (triangles and dashed lines). 



We begin our study of the connection by examining 
the cubic structure. The solid curve with filled circles 
(black online) in Fig. [s] is the inverse susceptibility calcu- 
lated for temperatures down to T = O.OSeV for the case 
U — 6eV, J — leV and carrier density n — 1.5 using 
the DOS shown in top panel of Fig. |4j which has the 
bandwidth W ~ 2.4eV. The inverse susceptibility curve 
is seen to deviate from the Curie form at low tempera- 
ture. The extrapolation to zero of the high temperature 
linear regime implies a Tc < 0. 

We attribute the flattening out of the susceptibility 
curve to the onset of Fermi-liquid coherence. To ver- 
ify this, and to gain additional insight into the na- 
ture of spin correlations in this model, we computed 

2 

the Wilson ratio Rw — (in our conventions 

3 7 

the dimensional factors — = 1). Here 7 is the co- 
Ms 

efficient of the linear specific heat 7 ^ 



dT 



T=0 



— Tr[i^i?Z~^].'^ We estimated the linear coefficient of 

the specific heat from the density of states vp at the 
Fermi level and calculated renormalization factor Z, 
which are obtained from the measured imaginary time 



Green's function via vp 



G{t = (i/2) and from 



the Matsubara frequency self energy via Z ^ sa 1 
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ImE(wi = 37r//3) - IinE(wo = t^IP) 



(not shown). 



The soHd curve with open circles (red onHne) in Fig. 
shows that, for the DOS in the top panel of Fig. |4]with 
U — 6cV and J — leV, the Wilson ratio extrapolates 
to the value Rw = 2 as temperature T — > 0. The 
value Rw = 2 is expected for a Kondo lattice with a 
low quasiparticle coherence scale but no intersite corre- 
lations, while a system with strong ferromagnetic cor- 
relations would be characterized by an Rw ^ 2. We 
therefore conclude that there is no evidence for ferro- 
magnetism in the cubic structure at C/ = 6, J = leV and 
W = 2.4eV and that the flattening of the x~^{T) curve 
indicates the onset of the Fermi liquid coherence. 

We now turn to the effects of the GdFeOs distortion. 
For simplicity of presentation, we focus mainly on the 
case 9 = (f). Fig. [4] shows the evolution of DOS with 
tilting angle. We see that, as the tilt angle is increased, 
the position in energy of the lowest density of states peak 
shifts down in energy. 
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FIG. 6: (Color online) Temperature dependence of inverse 
susceptibility computed for model Lai_a;Sra;V03 computed 
for tilt angles 6 = (fi from 8° to 18° computed at densities 
n = 1.55 (black circle lines) and n = 1.95 (diamond dashed 
lines) and interaction parameters U — 6eV, J = leV from 
Eq. |T]) with Hkin derived from MLWF fits to band struc- 
ture. The diamonds and circles are data points measured by 
DMFT. The lines are fitted from data points. Left column: 
plot over wide temperature range; right column: expanded 
view of smaU region. The vertical dashed line marks zero 
temperature. 

Fig. |6] shows the evolution of the inverse susceptibility 
with increasing amplitude of GdFeOa distortion at f7 = 
6eV and J — leV for two carrier density values n = 1.55 
and n = 1.95. For n = 1.55, ferromagnetic order can 



be observed starting a,t 6 = (j> = 14°. In contrast, for 
n = 1.95, the Curie temperature is nonzero within errors 
only for 6 = (p > 16° . 
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FIG. 7: (Color online) The magnetic phase diagram in the 
space of carrier density n (x-axis) and tilt and rotation an- 
gles 6 = (j) { 3/-axis) for the solid solution Lai-ajSr^^VOs with 
on-site interactions U — 6eV, J = leV. The white region indi- 
cates paramagnetic order (Tc < 0.004eV), the colored regions 
denote ferromagnetic order with the transition temperature 
Tc given by the scale bar (color on-fine). The solid curve is 
the 9 = 9{n) curve estimated for Lai-^SrajVOs as described 
in Appendix A. The estimated phase boundaries for U — 4 
(circle and dashed line, blue on-line) and 5eV (squares and 
solid line, green on-line) with J = leV are also plotted. 

The shaded areas (color online) in Fig. 7| show the 
phase diagram resulting from a detailed study of the de- 
pendence of Curie temperature on carrier density and 
tilt (rotation) angle at U = 6eV and J = leV obtained 
from curves such as those shown in Fig. [6] The expense 
of the computation and the uncertainties inherent in our 
extrapolation means that the phase boundary is not pre- 
cisely determined. We regard Tc < 0.004eV as consistent 
with Tc = within uncertainties. The width of the strip 
separating Tc = 0.004 and 0.008 gives a measure of the 
error bars on the locations of the phase boundary. The 
divergence in tilt angle required to obtain a nonzero tran- 
sition temperature as n — > 2 may be understood from the 
fact that for the correlation strength considered here the 
n = 2 compound is a Mott insulator. Ferromagnetism 
is favored by metallic motion of the carriers and is sup- 
pressed in proximity to the Mott insulating phase which 
for the parameters we consider is antiferromagnetically 
ordered. To summarize, we see that for these values of U 
and J obtaining a ferromagnetic state requires two con- 
ditions: (1) large hole doping away from LaVOa and (2) 
large GeFeOs-type distortion. 

We now consider the physics behind the dependence of 
transition temperature on tilt angle. Inspection of Fig. [4] 
shows that increasing the magnitude of the GdFeOs both 
reduces the bandwidth (from 2.4eV at = = to 
1.8eV at the critical angle) and moves the position of 
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the lowest density of states peak to lower energy. To in- 
vestigate the relative importance of the two effects we 
have performed computations for the cubic structure us- 
ing U = 6eV, J = leV but with reduced bandwidth 
W — 1.6 eV. The dashed curves in Fig. [sjshow the Wil- 
son ratio Rw (open triangles, red on-line) and the inverse 
susceptibility (filled triangles, black on-line) for the 
smaller bandwidth W = 1.6eV. The susceptibility curve 
indicates that the intercept is increased relative to the 
larger bandwidth case, but remains negative. The sharp 
downturn in the Wilson ratio at the lowest temperatures 
suggests that the decrease in intercept is a consequence 
of a lowered Fermi liquid coherence scale and does not 
indicate stronger ferromagnetic correlations. We there- 
fore conclude that for reasonable values of the bandwidth 
and correlation strength the crucial factor for ferromag- 
netism is the position of the density of states peak. For 
very large values of the Hund's coupling J(> 3eV), fer- 
romagnetism may occur over a much wider parameter 
range, as seen in previous Bethe lattice studies .1^^^^^ 




DOS maximum relative to the Fermi level (eV) 



FIG. 8: (Color online) The dependence of Curie temperature 
Tc on van Hove peak position with respect to the Fermi level 
of fully spin-polarized state calculated from Eqn. ([3| with 
t = 0.264eV, t' /t changes from —0.3 — >■ 0.3. On-site interac- 
tions U = 6eV and J = leV. The dashed line marks the zero 
temperature. 

The importance of the position of the density of states 
peak was previously noticed in the context of the one- 
band Hubbard model.'^HIll However, in that case ferro- 
magnetism was only found if the density of states peak 
is far below the Fermi level of the paramagnetic state. 
In the present case inspection of Fig. |4] shows that it is 
only necessary for the density of states peak to be not 
too far above the Fermi level. We believe that the differ- 
ence arises from the "double exchange" physics of Hund's 
coupling in partially filled d shells. The Hund's coupling 
favors high spin states, which means that hopping be- 
tween two sites is optimal if the spins are parallel and is 
suppressed if they are not parallel. This strongly favors 



ferromagnetism. A natural question is how far above the 
Fermi level can the density of states peak be and still 
support ferromagnetism. For a reasonable range of J 
(^ leV, see Fig. [s]) we find that a good rule of thumb is 
that ferromagnetism occurs if the density of states peak 
lies at or below the Fermi level of the fully polarized fer- 
romagnetic state. This answer is clearly not universal 
since model system studieJ^^^ti^ indicate that increasing 
J to very large values favor ferromagnetism even if the 
density of states peak lies very high in energy. 

We have also studied selected densities at the smaller 
correlation strengths U = 4eV and U — 5eV. Estimates 
of the resulting phase boundaries are shown as dashed 
and solid lines in Fig. [Tj We see that for intermediate n 
the phase boundary is only weakly dependent on U and 
n. For carrier concentrations near n = 2 weakening the 
correlations to move the system out of the Mott phase 
greatly increases the range in which ferromagnetism is 
observed. These calculations are the basis for the three 
schematic phase diagrams shown in Fig. [T] 



IV. THE EFFECT OF OXYGEN BANDS 

In early transition metal oxides, the bands that cross 
the Fermi level are well separated from other bands and 
are of mainly transition metal d character, arising from 
t2g symmetry orbitals. In many body calculations, it is 
common to focus only on these bands, which are treated 
as a multiba nd Hubbard model, while all other bands 
are neglected.^ ^ However, it is generally believed that 
the fundamentally correct model of the transition metal 
oxides should be based on atomic-like d orbitals coupled 
to weakly correlated oxygen p states.'^ In this section 
we examine the magnetic phase diagram of the vanadate 
system by applying the methods described in previous 
sections to the "pd" model which describes this situation 
and comparing the "pd" model results to those from the 
"d only" model. 

The pd model has two important energy scales: the 
on-site interaction U which as before is the energy cost 
for changing d occupancy and the charge transfer energy, 
in other words, the energy cost for an electron to hop 
from a ligand to a transition metal atom. The charge 
transfer energy may be defined in different ways but the 
correct value is not well established, either from ab ini- 
tio considerations or from experiment. The important 
points for our subsequent discussion are that the physics 
depends sensitively on the charge transfer energy as well 
as on U and that for reasonable U a metal-insulator tran- 
sition may be driven at integer band filling by varying the 
charge-transfer energy.l^ We will see that, as was found 
in the d only model, the ferromagnetic phase boundary 
depends on whether the parameters are such as to place 
LaVOs on the metallic or insulating sides of the metal- 
insulator phase diagram. 

The Hamiltonian describing the p and d states of a 
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FIG. 9: (Color online) Inverse susceptibility vs. tempera- 
ture for three different filling values 1.55, 1.90 and 2.00 for 
cubic (left column) and GdFeOs-distorted LaVOs with tilt 
angle 6 = tj) = 16° (right column). (a,b) pd model with 
large A (approximately that predicted by standarcP^ dou- 
ble counting correction); (c,d) pd model with small A; (e,f) 
d only model. On-site interactions U = 6eV, J = leV. For 
cubic case, Ajorge = 10.97eV, Asmaii = 2eV. For tilted case, 



Ajar 



lOeV, A, 



= OeV. The vertical dashed lines mark 



zero temperature. 



perovskite transition metal oxide has the fornP^l 



H — Hd + Hr, + H) 



hyb 5 



(6) 



where Ha = H^,„ + HonsUe with describing bare 

on-site energies and electron hoppings between d orbitals 
and Honsite as in Eq. ([4|; Hp describes the dispersion 
of the oxygen p bands in the absence of p hybridiza- 
tion; Hiiyii is the hybridization between d and p orbitals. 
We take Honsite to have the usual rotationally invariant 
Slater-Kanamori form with U = 6eV and J = leV. The 
band term Hkin = H^^^^ + Hp + Hhyi, has the same generic 
form as Eq. ^ but now the space of orbital indices a, /3 
is expanded to label both vanadium d and oxygen p or- 
bitals. The corresponding ^^^^^^(k) is generated using 
MLWF methods with a wide energy window including 
both d and p bands (see the MLWF procedure in Sec. 
and the A ppendix [b] ) . 

However j ^^l^^l^^l the Hartree shift arising from Honsite 
means that the d level energy obtained from the 
MLWF procedure must be renormalized by a "double 
counting correction" which we denote by A: = — A. 
We consider two values for A: one, which we refer to 
as ^large ^ lOcV, is essentially the value obtained by 
applying the standard^ double counting correction to 
basic band theory and is such that LaVOa is predicted 
to be metallic, in contradiction to experiment. The other 
value, which we refer to as A^maH ~ — 2eV, is such that 



the material is insulating at n = 2 in agreement with ex- 
periment. 

The resulting model is solved using single-site DMFT 
as described in Sec. |lTj but with one important addition. 
The full p d manifold includes V Cg orbitals. While the eg 
energy lies above the Fermi level, so that the e^-derived 
antibonding bands are empty, some of the filled bands 
are Cg-oxygen bonding states which have a small but non- 
zero eg content. One must therefore solve a 5-orbital im- 
purity model. Performing this calculation in full general- 
ity would be prohibitively expensive. We therefore follow 
standard procedure and treat the eg orbital contribution 
to the impurity model in a Hartree approximation. 

Fig.[9]shows representative computations of the inverse 
susceptibility for cubic and tilted systems at several car- 
rier concentrations. The top row (panels a and b) dis- 
plays pd model results obtained for the standard dou- 
ble counting correction (so that LaVOa is wrongly pre- 
dicted to be a metal) while the middle rows (panels c 
and d) show results obtained if the double counting cor- 
rection is tuned so that the calculation places LaVOa in 
the Mott/charge-transfer insulating regime of the phase 
diagram. The bottom two panels (e and f ) present d only 
model results for comparison. In the d only model the U 
value is such as to place the n — 2 (LaVOa) material in 
the Mott insulating region of the phase diagram). The 
transition temperature estimates obtained by linearly ex- 
trapolating the curves to are given in Table |l] 
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0.1167 
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cubic, d only model, U — 


4eV 


-0.2347 


-0.1368 


-0.1250 


tilted, d only model, U — 


4eV 


0.0267 


0.0345 


0.0366 



TABLE I: Values for Curie temperature Tc (in eV) for each 
case considered in Fig. [9] together with results for U = 4eV 
for three different fillings n = 1.55, 1.90 and 2.00. All compu- 
tations are for J = leV; except where indicated, U — 6eV is 
used. 

Examination of the results in Table. U shows that the 
qualitative trends are the same in the pd and d only model 
calculations. In particular, in both models increasing the 
tilt angle increases the tendency towards ferromagnetism. 
However, significant differences are visible; in particular 
the pd model has a significantly greater tendency to fer- 
romagnetism than does the d only model and (especially 
in the small-A case) the differences are more pronounced 
for the cubic than for the tilted structure. 

We believe that there are two main origins for the dif- 
ferences. First, in the small A cubic system case, the 
change in the charge transfer energy relative to band the- 
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ory affects the density of states, moving the peaks closer 
to the Fermi level while for the tilted structure the shift 
in charge transfer energy does not change the peak po- 
sitions as much (see Appendix C). Second, and perhaps 
more important, the Cg state occupancy arising from the 
p — d bonding bands (omitted in the d only model) in- 
creases the effective moment on the d site, thereby en- 
hancing the tendency towards magnetism. This effect is 
more pronounced in the larger A (smaller charge transfer 
energy) case, because the p d mixing is larger. 
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FIG. 10: (Color online) Inverse susceptibility vs. temperature 
for cubic (blue color) and 16°-tilted (red color) structures at 
filling n = 2.0. (a) Comparison for pd model with (solid lines) 
and without (dashed lines) eg band spin polarization, (b) 
Comparison between d only model (solid lines) and pd model 
without Eg band spin polarization (dashed lines) . Parameters 
for calculations with pd model are the same as in Fig. [9^,b. 
Calculations with d only model use U = 4eV < Uc,J~ leV 
for both structures. 



To demonstrate this point we present in Fig. 10 i,b cal 



culations of the inverse susceptibility under different con- 
ditions. The curves in the left panel compare calculations 
in which the Cg occupancy is frozen at the spin unpolar- 
ized values (dashed lines) and calculations in which the eg 
orbitals are treated within the Hartree-Fock approxima- 
tion as described above. We see that the feedback from 
the polarization of the Cg orbitals makes a significant con- 
tribution to the transition temperature. The right pan- 
els show that calculations performed in the "frozen e^" 
model agree reasonably well with the corresponding cal- 
culations in the d only model for U = AeV < Uc- Further 
analysis is given in Appendix D. 



V. CONCLUSIONS 

In this paper we have investigated the conditions under 
which ferromagnetism might be observed in bulk solid 
solutions derived from LaVOa. Our theoretical studies 
used realistic band structures derived from maximally- 
localized Wannier function fits to density functional (typ- 
ically PWscf/GGA) band structure calculations, along 
with the single-site DMFT approximation which contains 
the local physics of partially filled orbitally degenerate 
d orbitals with strong Hund's coupling. Our main the- 



oretical finding is that for physically reasonable values 
of the Hund's coupling ferromagnetism is driven mainly 
by density of states effects, being favored by a density of 
states peak lying not too high in energy. The dependence 
on correlation strength and carrier concentration is rela- 
tively weak except that if the correlation strength is large 
enough to drive a Mott transition to an antiferromag- 
netic Mott/charge transfer insulator, ferromagnetism is 
suppressed in the vicinity of the insulating state. 

An association between ferromagnetism and the den- 
sity of states was previo usly no ticed in studies of the 
one band Hubbard modeLfSEHlll However, there is a sig- 
nificant difference between ferromagnetism in the model 
studied here and that found in the single-band Hubbard 
model. In the single-band Hubbard model, ferromag- 
netism only occurs when the DOS peaks are located far 
below the chemical potential, and indeed very close to the 
lower band edge^ In the orbitally degenerate, Hund's 
coupled systems studied here, the requirements on the 
position of the density of states peak are substantially re- 
laxed. We find that for physically reasonable correlation 
strengths, ferromagnetism can occur as long as the low- 
est density of states peak is at or below the Fermi level of 
the fully polarized ferromagnetic state. The importance 
of orbital degeneracy and Hund's coupling has been pre- 
viously noted in studies of models with a semicircular 
(Bethe-lattice) density of st&te^^^^ and in the periodic 
Anderson modeP^but the important role played (for rea- 
sonable values of J) by the near-Fermi-surface density of 
states peaks seems not to have previously been noticed. 

In the vanadate systems, the density of states is con- 
trolled by the amplitude of the GdFeOa distortion away 
from the ideal cubic perovskite structure. We find 
that increasing the distortion amplitude favors ferromag- 
netism. Indeed, if the U interaction were set to the un- 
physically small value U — 4eV so that LaVOa were not 
a Mott insulator, then the observed GdFeOa distortion 
would be large enough to put the material at or very near 
to the ferromagnetic phase boundary. One may then ask 
why doping does not induce ferromagnetism in bulk solid 
solutions such as Lai_a:Sra;V03. We believe that the an- 
swer is that in the physical system, the distortion ampli- 
tude and the carrier concentrations are both determined 
by the La/Sr ratio in such a way that, as x is varied, the 
distortion amplitude decreases so that the system follows 
a path in parameter space which remains outside the fer- 
romagnetic region (solid line in Fig. [?]). 

We remark that while the general experience has been 
that the single-site DMFT approximation provides a 
good qualitative representation of the physics, i n partic- 
ular of trends as material parameters are varied,El22l22li£l 
its quantitative accuracy in producing magnetic phase 
boundaries has not been established. It is possible 
that single-site DMFT underestimates magnetic order. 
The reports of ferromagnetism in vanadate superlattice^ 
where both the doping and the tilt angle may be locally 
large, suggests that the true ferromagnetic phase bound- 
ary may be shifted to lower tilt angle but still above the 
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bulk LaVOs 9{n) curve. Increases in computer power 
and improvements in algorithms may soon make cluster 
dynamical mean field studies of realistic systems with or- 
bital degeneracy and Hund's coupling feasible, and it will 
be important to use these methods to assess the reliabil- 
ity of single-site DMFT predictions. 

The second point of materials theory is the issue 
of what is the relevant theoretical model. Consistent 
with most th eoretical literature on early transition metal 
oxides J^^^Bs! focused mainly on a c? only "multiband 
Hubbard model" approach where the electrons in the 
near Fermi surface d derived bands were correlated by lo- 
cal U, J interactions. The results are in reaso nable accord 
with experiment. It has been arguecp2l2Sl2ll that a more 
generally valid description may be obtained from a model 
in which correlations are applied to atomic-like d orbitals 
which are hybridized to the other orbitals, in particular 
to the oxygen p orbitals. In this "pd" model approach, a 
key parameter is the renormalized d level energy, which is 
shifted from the band theory value by a "double counting 
correction" for which there is no generally agreed value. 
If the standard expressioiP^ for the double counting cor- 
rection is used, in the single-site DMFT approximation 
LaVOa is found to be a ferromagnetic metal rather than 
an antifcrromagnetic insulator and indeed the solid so- 
lution Lai_a;Srj.V03 is predicted to be ferromagnetic for 
a wide range of x, in reasonable correspondence to the 
small U calculation in the d only model. If the double 
counting contribution is shifted such that the model for 
LaVOa is in its Mott insulating regime, then the result- 
ing phase diagram is very similar to that shown in Fig. [7] 
for the d only model. Quantitative differences arise from 
the changes in density of states arising from the large 
renormalization of required to make the n = 2 case a 
Mott insulator and from the contribution of the (small 
but non-negligible) eg occupancy arising from the the p-d 
bonding states. Determining the appropriate theoretical 
approach for the early transition metal oxide is an im- 
portant open question. 

Our findings suggest several routes to creating ferro- 
magnetism in artificially fabricated systems. The gen- 
eral aim should be to manipulate the band structure so 
as to move density of states peaks close to the Fermi 
level. Inducing octahedral tilts by compressive strain is 
an important route. Further, in early transition metal 
oxides, increasing the p d hybridization clearly increases 
the effective d moment and this provides a self-consistent 
amplification of the Curie temperature. Therefore ma- 
nipulation of the p d energy difference can be used to 
control Tc. 

We also observe that some transition metal oxides such 
as the Ru-based perovskites and Ruddlesden Popper ma- 
terials Sr„4.iRu„03„_|_i involve holes in the bands and 
in a qualitative sense may be thought of as the particle- 
hole transforms of the models studied here. In this case 
for the physical sign of t' the van Hove peaks are on the 
other side of the Fermi level, suggesting that the theoret- 
ical models will be more likely to exhibit ferromagnetism. 



However, increasing the tilt angles shifts the peaks in the 
wrong direction, reducing the tendency to magnetism. A 
more detailed investigation of this physics is in progress. 
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Appendix A: 6{x) Curve for Lai a^Sr^jVOa 

The relation between tilt angle 9 and the doping level 
X is obtained based on the fact that the (pseudo) cubic 
lattice constant of Lai-j^Sr^-VOa is linearly dependent 
on i.Eal Bulk SrVOa has a^^° = 3.84A,l2Sl bulk LaVOg 
has ttp^^ = 3.92A,[22l hence for doped material, ap{x) = 
af^Ox + a^^O(l-x). 

We further assume that the effective V-0 distance d 

also has similar relation, d(x) — dsvo^ + c?Lyo(l ~ x)- 

^svo 

Bulk SrVOa is cubic, thus dsvo = = 1.92A; while 

dLVO = 2.00A.f2Sl The angle 6 is calculated from Up = 
2c? cos 0, and with n = 2 — a; we obtain the 9{n) curve 
shown in Fig. [7] 

Appendix B: Band Calculations and Wannier Fitting 

Our procedure for obtaining the band theoretic contri- 
bution to the local Green's functions is as follows. First, 
given a rotation angle (j) and tilt angle 9 we construct 
the 4-unit-cell structure corresponding to a GdFeOs-type 
distortion of the cubic p erovs kite. The PWscf code of 
QUANTUM ESPRESSCEnmi is then used to calculate 
the band structure for that system. Finally, we use 
the Wannier90 implementatiorP^l of the MLWF proce- 
dure with an appropriate energy window to generate our 

For the band structure calculations with PWscf, we 
used a cutoff energy Ecutoff = 30Ry w 408eV and a 
10 X 10 X 10 Monkhorst-Pack k-mesh. The MLWF proce- 
dure involves the choice of both an overall energy window 
and a "frozen" window within which the MLWF bands 
are forced to coincide with the DFT bands. For MLWF 
fitting with Wannier90, we set the overall energy win- 
dow for the d bands to run from — 1.5eV to 6eV with 
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respect to the chemical potential. For larger octahedral 
rotations, the Cg and t2g bands are well separated and 
choosing the frozen energy window to run from — 1.5eV 
to 1.5eV represents the t2g bands well. For smaller ro- 
tations, band overlap becomes important and the upper 
cutoff of the frozen energy window must be reduced to 
avoid overlap with the eg bands; this means that the 
high-lying (unoccupied) part of the t2g bands is not per- 
fectly represented, but these issues do not affect our main 
results, which concern the location of the ferromagnetic 
phase. If oxygen bands are included in the calculation, 
the energy window is enlarged from —10 to 6eV, while the 
range —10 — ?> 1.5eV is set as the frozen energy window. 



Appendix C: Charge Transfer Model 

In the charge transfer {pd) model, the value of the 
charge transfer energy A affects the ferromagnetic phase 
boundary and also changes the band structure. In Fig.[TT] 
we demonstrate the effect of varying A on the band struc- 
ture of cubic LaVOa . 

We use the same parameters as in Sec. |IV| for the 
cubic structure and employ a Hartree-Fock approxima- 
tion to derive the spin/orbital unpolarized DOS. Fig. 11 
shows such DOS for A = 2eV and A = 10.97eV. Within 
Hartree-Fock calculation, oxygen p bands are located at 
the same positions as in corresponding DMFT results, 
so the p d hybridization may be expected to be simi- 
lar in the two calculations. We therefore consider the 
Hartree-Fock DOS as the "non-interacting" DOS for the 
pd model, including charge transfer effects. 
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FIG. 11: (Color online) Hartree-Fock orbital/spin unpolarized 
t2g band density of states with U = 6eV and J = leV for two 
cases of cubic structure used in Sec. |IV| /^smaii = 2eV and 
Ajarge = 10.97eV. The dashed vertical line marks the Fermi 
level. 



The two curves shown in Fig. 11 demonstrate that 
changing A changes both the position of the van Hove 
peaks relative to the Fermi level and the bandwidth. 
For Aiarge = 10.97eV, the bandwidth is = 2.58eV 



and the ratio between next nearest neighbor and nearest 
neighbor hoppings is t' /t = 0.23. For Agmaii = 2eV, 
the bandwidth W of the antibonding band is smaller 
(W — 2.03eV) and the van Hove peak moves closer to 
the Fermi level {t' /t — 0.18). These two factors act to 
make the model more ferromagnetic. 



Appendix D: Contribution of Cg Bands 

In present-day applications of the DFT+DMFT 
method to transition metal oxides, atomic-like d orbitals 
are defined via a projector or Wannier construction and 
are coupled to other bands (mainly oxygen p, in practice) . 
The resulting theory differs from the multiband Hubbard 
models often used^ in two ways. First, the explicit inclu- 
sion of oxygen orbitals means that charge transfer physics 
in the sense of Ref . E51 mav play a role. Second, the full d 
manifold is included, in contrast to the multiband Hub- 
bard approach where only the frontier orbitals (in the 
present case, the states) are treated. This latter differ- 
ence is important for the vanadates because even though 
the Cg orbitals are notionally empty, hybridization with 
the oxygen bands leads to a non- vanishing density of Cg 
electrons. 

A fully interacting treatment of the entire d manifold is 
too expensive to be computationally practicable. In this 
paper we take the virtual occupation of the Cg orbitals 
into account via a Hartree approximation. In each iter- 
ation of the DMFT self consistent calculation, eg occu- 
pancies obtained from the previous iteration are used to 
calculate the eg feedback to the t2g bands via the Hartree 
approximation to the inter-orbital self energy: 



J2 [(U - 2 J)n,s + {U- 3 J)n,,] . (Dl) 



bands 



The two panels of Fig. [TO] show that it plays an impor- 
tant role in the magnetic phase diagram. In what follows 
we give analytic Stoner-style arguments explaining these 
results. 

To estimate the contribution of eg bands to the polar- 
ization, we first observe that although the eg orbitals are 
notionally empty, hybridization with the oxygen states 
means that the eg occupancy is non-vanishing; roughly 
if the p Bg energy difference is Ag and the p — d hopping 
is tpd then 



a:- 



(D2) 



A perturbation to the eg energy of magnitude h thus 
leads to a change in occupancy 



J tpdh 

<^g \2 ' 

e 

SO that we expect a susceptibility given by 

a!' 



Xe 



(D3) 



(D4) 
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The absence of eg character at the Fermi level and the 
fact that the eg — p energy difference (Ae) is larger than 
that of the t2g — p energy difference suggests that Xeg 
will be smaller than the band susceptibility, but not by 
orders of magnitude. 

We will assume that the intra-i2g contributions lead to 
a Xt2g = ^/(T — Tq) which would diverge at some tem- 
perature To . In the presence of a magnetic field this Xt2g 
leads to a polarization of the t2g orbitals which, via the J 
interaction, affects the Cg-oxygen p energy difference and 
thereby leads to a polarization of the eg orbitals which 
in turn feeds back on the t2g states. Mathematically, we 
distinguish the eg and t2g orbitals and write me^ tj^ as 
the magnetization of one orbital 

mt^^ ^ Xt2gih + 2.hne^), ^^^^ 



Solving these equations gives 



^9 ~ ~ {XegY'-J (D6) 

^A-\T^To)~&J\,g. 



Thus the effect of the eg component of the bond- 
ing bands is to shift the Curie temperature by AT « 
(yj'^Axeg- From Fig. [lOj we estimate A « 3 so even 
though although Xeg isTikely to be small, the overall ef- 
fect may be non-negligble. 
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